Multiscale finite element and domain decomposition methods for high-contrast problems using local spectral basis functions
نویسندگان
چکیده
In this paper, we study multiscale finite element methods (MsFEMs) using basis functions that provide an optimal convergence in domain decomposition methods. We consider second order elliptic equations with highly variable coefficients. Both multiscale finite element and domain decomposition methods (considered here) use coarse spaces to achieve efficiency and robustness. In MsFEMs, the spatial variability of the media affects the convergence rate (see e.g., [14]). In domain decomposition methods, high variability of the coefficients within coarse regions affects the number of iterations required for the convergence. In this paper, we use coarse spaces designed for high-contrast problems that provide an optimal convergence in domain decomposition methods. These spaces are constructed locally. In particular, basis functions are constructed using solutions of a local spectral problem. In our previous work [18], we show that using these coarse spaces one can construct preconditioners such that the condition number of the preconditioned system is independent of the contrast. In this paper, these coarse spaces are used in MsFEMs to solve elliptic equations with high-contrast heterogeneous coefficients on a coarse grid. Our numerical results show that MsFEMs with coarse spaces constructed via local spectral problems are more accurate compared to multiscale methods that employ traditional multiscale spaces, e.g., with linear boundary conditions. However, we argue that these coarse spaces based on local eigenvalue problems are not sufficient to capture the fine-scale behavior of the solution very accurately. In particular, these coarse spaces identify high conducting regions accurately, but they may not capture the detailed behavior of the solution in other regions, such as between high-conductivity regions. The latter can be important for the ∗Department of Mathematics, Texas A & M University, College Station, TX 77843 ([email protected]) †Department of Mathematics, Texas A & M University, College Station, TX 77843 ([email protected]) ‡ExxonMobil Upstream Research Company, P.O. Box 2189, Houston, TX 77252 ([email protected])
منابع مشابه
Domain decomposition preconditioners for multiscale problems
In this paper, we study domain decomposition preconditioners for multiscale elliptic problems in high contrast media. We construct preconditioners such that the condition number of the preconditioned system is independent of media contrast. For this purpose, multiscale spaces for the interpolation on the coarse grid is developed using a local weighted spectral problem. A main observation is tha...
متن کاملA Multiscale HDG Method for Second Order Elliptic Equations. Part I. Polynomial and Homogenization-Based Multiscale Spaces
We introduce a finite element method for numerical upscaling of second order elliptic equations with highly heterogeneous coefficients. The method is based on a mixed formulation of the problem and the concepts of the domain decomposition and the hybrid discontinuous Galerkin methods. The method utilizes three different scales: (1) the scale of the partition of the domain of the problem, (2) th...
متن کاملGeneralized Multiscale Finite Element Method for Elasticity Equations
In this paper, we discuss the application of Generalized Multiscale Finite Element Method (GMsFEM) to elasticity equation in heterogeneous media. Our applications are motivated by elastic wave propagation in subsurface where the subsurface properties can be highly heterogeneous and have high contrast. We present the construction of main ingredients for GMsFEM such as the snapshot space and o in...
متن کاملA Systematic Coarse-Scale Model Reduction Technique for Parameter-Dependent Flows in Highly Heterogeneous Media and Its Applications
In this paper, we propose a multiscale approach for solving the parameter-dependent elliptic equation with highly heterogeneous coefficients. In particular, we assume that the coefficients have both small scales and high contrast (where the high contrast refers to the large variations in the coefficients). The main idea of our approach is to construct local basis functions that encode the local...
متن کاملImplementation of a Mortar Mixed Finite Element Method using a Multiscale Flux Basis
This paper provides a new implementation of a multiscale mortar mixed finite element method for second order elliptic problems. The algorithm uses non-overlapping domain decomposition to reformulate a fine scale problem as a coarse scale mortar interface problem, which is then solved using an iterative method. The original implementation by Arbogast, Pencheva, Wheeler, and Yotov, Multiscale Mod...
متن کامل